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I. INTRODUCTION 

The Hamiltonian iL of a physical system represents its energy^which is a real observable. It is therefore required 
that the expectation values of the quantum operator H be real [J. This can be guaranteed by imposing that H be 
Hermitian, H = H\ as it is known that the spectrum of a Hermitian operator is real and its eigenvectors form a 
complete orthogonal set 0 . 

It is on the other hand known that Hermiticity is not a necessary condition for a real spectrum Q : a large number 
of one-dimensional non-Hermitian potentials, both real and complex, invariant under the simultaneous actions of the 
parity P (space reflection) and time reflection T operators have been found to admit energies that are real and 
discrete. 

The matter is not a idle one, as non-Hermitian PT-invariant operators find applications in many areas of theoretical 
physics: “optical” or “average” potentials in nuclear physics quantum field theories scattering problems Q, 
localization-delocalization transitions in superconductors Q , detraction of atoms by standing light waves 0 , as well 
as the study of solitons on a complex Toda lattice 0 - 

Unfortunately, PT-invariance is neither necessary nor sufficient to ensure the reality of the spectrum; however, 
it has been conjectured Q that PT invariant Hamiltonians possess real discrete eigenvalues if the PT symmetry is 
unbroken i.e. if the energy eigenstates are also eigenstates of the operator PT. When the PT-symmetry is broken and 
the Hamiltonian is real there instead are energy eigenvalues that are complex conjugate pairs. However, no general 
condition has been found for the breakdown of the PT-symmetry. 

In this contest, it has been pointed out that a necessary, but not sufficient, condition for the spectrum to be real and 
discrete is the 77 -pseudo-Hermiticity, rjHrj~^ = , of the Hamiltonian, where 77 is a Hermitian linear automorphism 

0 . The property is also known as selfadjointness in an indefinite inner product space, see [IMl. The eigenvectors 
of H are in this case 77 -orthogonal, i.e. they are orthogonal according to the 77 -distorted inner-product < >■ 

Several PT-symmetric potentials have been found to be P-pseudo-Hermitian [l^ and classes of non-Hermitian 
Hamiltonians -both PT-symmetric and non-PT-symmetric- appear to be pseudo-Hermitian under rj = where 
0 G R and P = , (^=l)is the momentum operator (the transformation generated by 77 is an imaginary shift: 

r}xrj~^ = X -\- iO. npri~^ = p) [l^ . or 77 = where t{x) is a function of x (the transformation is a complex 

gauge-like one) [I 4 . 

It is thus still a case by case procedure to check whether the eigenvalues of an operator are all real. This does not 
usually cause big practical problems when dealing with a single operator. The situation changes when we have to 
consider classes of operators. Procedures have been developed for families of operators acting on spaces with finite 
bases, see e.g. Ref. [Tsi l whose author considers a one parameter family of PT-symmetric matrices M{e), with a 
perturbation parameter e G R which destroys Hermiticity while it respects PT-invariance. 

Here we consider the case when for numerical simulations an operator H acting on a space with an infinite basis 
needs to be truncated and study the rate of convergence to their asymptotic value of those eigenvalues that for 
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truncated matrices may happen to be non-real. The operator H = i?[o,oo) is given by a non-symmetric Jacobi matrix 


H, 


(oo —i>o 

bo ai bi 

bl 02 62 


[0,oo) 


V 


&2 03 ■ ■ 


(LI) 


with bounded, real sequences {aj)^Q, {bj)^Q, the sequence {bj)^Q being additionally strictly positive. Its finite 
truncations are of the form 


^[0,n] 


/Oo -bo 
bo oi bl 

bl 02 




1 

1 O^ j 


( 1 . 2 ) 


and 


?7 = diag (-1,1,1,...). 

Due to the fundamental theorem of Pontryagin [T^ each the operators LJ[o,n] ha-s, generically, either a unique single 
eigenvalue A„ on the real axis with the eigenvector /„ satisfying {fn\vfn) < 0 or a single couple of complex conjugate 
eigenvalues A„ £ C+, A„ £ C“ (to avoid confusion with the conventions used in some of the papers we quote, we note 
that here and in the following {x\y) always denotes the usual inner product -either in C” or in P'- linear with respect 
to the second variable). The remaining part of the spectrum of LJ[o,n] is real. The same is true for the spectrum of the 
infinite matrix iL[o,oo] with the eigenvalue Aoo, see Section [H] for details. The character of the convergence A„ —>■ Aoo 
is the main topic of our paper. 

Our approach makes use of analytic representations of the function 

^[0,oo)('^) {^01 (.^[ 0 , 00 ) ; (i-’ii) 

which contains the full information about the spectrum of LJ[o,oo) of its {n — 1/n] Fade approximants 


W[o,„](z) = - (eo|(iL[o.n] - ^) ^eo). 


(1.4) 


In particular, A„ (Aoo) is a pole of 'm[o,n\{z) (™[o,oo]; respectively) and it can be characterized in analytic terms. Due to 
the locally uniform convergence of m[o,n\ to TO[o^oo) the sequence (A„)))Tq converges to Aoo ('Corollary III. 51) . Our 
main interest is the rate of this convergence. In particular we show its dependence of the placement of the eigenvalue 
Aoo in C’*' U K. 

Our paper is organized as follows: 

• We give various analytic representations of the function dLU, choosing in particular as our starting point 


-1 

W[0.oo)(^) 


— flo — z J- 6 q 



dy{t) 
t — z’ 


where /i is some probability measure, see Theorem lII.il 


• In the case Aoo 4- [^ 81 ^ 4 ] show that the convergence rate of A„ to Aoo is exponential, with the base of the 
exponent increasing with the distance of Aoo from [^ 3 ,^ 4 ]: see Theorem IIII.3I below. 


• However, the A„’s tend to arrange themselves in branches spiraling into Aoo and some of these branches can get 
trapped in the real axis for a number of iterations no which can be relatively large when Aoo is close to [^ 3 ,^ 4 ]. 
We show examples with different numbers of branches and compute an estimate for no in Theorem lA.ll 


• In the case when Aoo £ [^3 j f 4 ] we build an example to show that the convergence rate is in general worse than 
exponential. 


• In the concluding remarks we review the possible cases from the numerical point of view. 















3 


II. HOLOMORPHIC REPRESENTATIONS OF THE m-FUNCTION 

We start with reviewing the spectral properties of the matrices i?[o,n] and i?[o,oo]- The matrix i?[o,n] is selfadjoint 
in the indefinite inner-product space with the fundamental symmetry given by rjn = [—1] ®In (? 7 ra-pseudo-Hermitian). 
Consequently one of the following four possibilities applies: 

(i) is similar to a diagonal matrix with real entries, except two complex conjugate entries A„ £ C^, A„ £ C“. 

The eigenvectors /„,g„ corresponding to the eigenvalues A„,A„ of satisfy (/„|77„/n) = (5nbn5n) = 0, 

(/nl^nffn) 7^ 0. 

(ii) H[Q,n] is similar to a diagonal matrix with real entries and there is precisely one eigenvalue A^ with the corre¬ 
sponding eigenvector /„ satisfying (/„|77„/„) < 0. 

(iii) is similar to a block-diagonal matrix with all the blocks real and one-dimensional, except one block of the 
form 


(K 1 

'vO A„ 


with A„ £ M. 


The eigenvector /„ corresponding to the eigenvalue A„ of H[o,n] satisfies {fn\r]nfn) = 0. 

(iv) -ff[o,n] is similar to a block-diagonal matrix with all the blocks real and one-dimensional, except one block of the 
form 

/A„ 1 0\ 

I 0 An 1 I with An £ R- 
V 0 0 A„y 


The eigenvector /„ corresponding to the eigenvalue An of ^^[o,n] satisfies (fn\Vnfn) = 0. 

The cases (iii) and (iv) are non-generic, i.e. the set of all matrices ^^[o,n] for which one of them applies has measure 
zero. We refer the reader to [l^ for the full canonical form of matrices selfadjoint in indefinite inner-product spaces, 
which gives also a full description of the eigenvectors. We observe that the matrix i?[o,n] may jump back and forth 
with n among the four types above. 

The spectral properties of the infinite matrix i?[o,oo)i understood as an operator on are more tricky: we refer the 
reader to for a full description and for canonical models. Here we note only that again there are essentially 

two possibilities: 


(i’) i?[o,oo) is similar to an orthogonal sum of a bounded selfadjoint operator in a Hilbert space and a diagonal 
matrix with two complex conjugate entries A„ £ C"*", A„ £ C“. The eigenvectors fn,gn corresponding to the 
eigenvalues A„, A„ of i?[0,oo) satisfy (/n|??n/n) = (5n|?7n5n) = 0, {fnlVngn) 7^ 0. 

(ii’) The spectrum of if[o,oo) is real and H^[o,oo) has a (unique) real eigenvalue with the corresponding eigenvector /oo 
satisfying (/oo|?7/oo) < 0. 


In the (ii’) case the Jordan chain correspondin g to Aoo is again of length not greater than three. 

Now we specify the theory developed in [^, to the case we are dealing with in the present work. Besides 

the matrices iJ[o,oo) and H[o,n] defined in (II.II) and (11.21) . we shall use the following truncations of the matrix i?[o,oo) 




/ ai h 

bi ai 

V 


1 

^n—1 ) 


n = 1,2, 


(II.l) 


Furthermore, JJ[i,oo) '"’iU stand for the infinite, symmetric Jacobi matrix with (uj)^]^ on the main and on the 

second diagonals. Similarly to m and (ILI we define the functions 

m[i,„](2:) = (ei|(i?[i,„] - z)~^ei) , m[i,oo)(z) = (ei|(FJ[i,oo) - ■ (11.2) 


Here ej stands for the j-th vector of the canonical basis of We call the functions appearing in (|L3|, dLll) and 
(III. 21) the m-functions of the corresponding Jacobi matrix. We refer the reader to for a treatment of m-functions 
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of symmetric Jacobi matrices appearing in (III. 21) . The fuirctions m^i n] (n £ Z_|_) and TO[i_oo) £^re analytic in the open 
upper half-plane C'*'. The function m[o,oo) (^[o,™]) is analytic in the upper half-plane, except Aoo (A„, respectively). 
Moreover, the Schur complement argument provides the following crucial relations [^, 




1 

z- ao-bl TO[ 1 _„]( 2 )’ 


z £ C"*" \ {An}, n £ Z-)-, 


(11.3) 


™[0.oo)(2:) = 


z-ao-b'^ m[i_oo)(^)’ 


z£C+\{Aoo}. 


(II.4) 


Let us now recall the definition of the class A/i. By Afi we define the set of generalized Nevanlinna fuirctions with one 
negative square, that is the functions of one of the three forms 


{z-P)[z-pf^ 

(11.5) 

- ~T{z)^ 

( 11 . 6 ) 

(z - a)(z -a)(p{z}, 

(11.7) 

Nevanllina function, i.e. 

(f is holomorphic in C+ and maps C+ into 


subclass of A/i functions to be investigated in the present work: 

Theorem II.1 Letm be a meromorphic function in the open upper half plane. The following conditions are equivalent. 

(i) There exist Aoo £ U M, d £ R and a nontrivial Borel measure a having all moments finite and supported on 
an interval [^ 1 ,^ 2 ] such that 


m{z) 


1 

( 2 : Aoo)(.2: Aoo) 


(^z + d + 




(II. 8 ) 


(ii) There exist oq £ R, 60 > 0 and a nontrivial Borel probability measure p having all moments finite and supported 
on an interval [^ 3 ,^ 4 ] such that 


-1 _ ^,2 dp{t) 

— Oq — z + b^ I ^ ^ 


MZ) Jt3 

(Hi) There exist a matrix id[o,oo) of the form (II.II) with bounded entries Oj £ R, bj > 0, j £ Z_|_ sueh that 

m{z) = TO[o,oo)( 2 ) := - (eo|(id[o,oo) - z)~^eo} . 


(II.9) 


(II.IO) 


Furthermore, the parameters Aoo and d and the measure a in (i), ao,&o ond p in (ii), and aj,bj, j £ Z+ in (Hi) are 
uniquely determined; the numbers qq and bo in statements (ii) and (Hi) coincide and Aoo from statement (i) is the 
(unique) eigenvalue of nonpositive type of the operator FIfrom statement (Hi). 

Proof. The equivalence (ii)<t^(iii) is a consequence of equation (III. 41) and the classical theory which sets a corre¬ 
spondence between the functions mp and the Jacobi matrices iJ[i_oo)) see e.g. [ 2 ^,[ 2 ^. 

(iii)=J>(i) Let m = m[o,oo). From the construction in [2lj it follows that m belongs to the class Ni and hence it has 
one of the forms (IIL5I) - (IIL7I) . 

Furthermore, expanding the resolvent into a geometric series at infinity one sees that m necessarily possesses an 
asymptotic expansion at infinity 

m{z) = - (eo|(iL[o,oo) - z)-^eo) = ^ ^--, (H.ll) 
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with Sj G M (j = 1, 2^. ) (see [13, Ull)- Comparing the forms pi.Sp - dll.TI) with pi.lll) . one gets by the Hamburger- 
Nevanlinna theorem [Ig] that the function (p in (III.5I) - (III.7I) can be represented in the form 


p(z)=z+d+ r ( 11 . 12 ) 

7ti t-z 

where d G M, and cr is a measure with all moments finite. Furthermore, comparing the expansions of (III. 51) . (III. 61) and 
(III. 71) with (III. Ill) we can see that case pi.6l) applies. In consequence, 


m{z) 


1 

Aoo)('2^ Aqo) 


(^z + d + 




(11.13) 


Observe that the measure cr cannot be a finitely supported (trivial) measure, since bj > 0 for all j G Z_|_ and in 
consequence neither m[i^oo) nor m[o,oo) are rational functions. The uniqueness of the parameters Aoo and d and of 
the measure a follows from the theory of Afi functions, see e.g. [ 13 ]. The fact that Aoo is the nnique eigenvalue of 
nonpositive type of i7[o,oo) follows e.g. from Ref. [ 13 ] or [^. 

(i)^(ii) Using the algorithm proposed in [ 13 ] (see also [ 13 |), one can find that m defined by (111.131) can be uniquely 
represented as 


m(z) 


1 

z — oo — 5 q rni{z) ’ 


where toi(z) = 


is a Nevanlinna function with finite moments. 


(11.14) 

□ 


Remark II. 2 Already at this point we can say something about the influence of the the measure p (spectrum of the 
matrix i7[i_oo)) on the position of Aoo- 

1) Conditions (i) and (ii) in Ref. [HI tell us that Aoo G [fa, < 4 ] if and only if 


[ |<-Aoo| ‘^dn{t)<bQ^, ao-Aoo+&o/ (^-Aoo) ^dp{t) = 0. (11.15) 

Jt3 Jt3 

(In particular, since bo is strictly positive, for the first of these conditions to be true, \t — Aoo|~^ needs to be a /i- 
integrable function; the second condition is just the specialization of eq. (111.91) to the case 2 : = Aoo and we introduce 
it here to fully characterize Aoo itself). It follows that if the measure p is sufficiently dense Aoo cannot be on [^ 3 ,^ 4 ], 
i.e. jjL “repels” the point Aoo- This is the case for Examples 1111.41 and 1111.51 below. 

2) If instead p has gaps, these gaps tend to trap Aoo- To see this, let’s assume that the spectrum of -ff[i,oo) has 
a gap {to,to) C \tz,t4\ and that uq G {to,to). From the definition of i7[o,oo), eq. (11.11) . it’s obvious that for bo = 0 
the point oq is an eigenvalue of H^o.oo) with the corresponding eigenvector / satisfying {f\r]f) < 0 , i.e. oq = Aoo- 
We now increase bo = 0; applying Rouche’s theorem to — l/m[o,oo) and remembering that if Aoo ^ K then Aoo is also 
an eigenvalue, we see that Aoo moves along the real axis until it meets anoth er p art of the spectrum of i7[o,oo) (that 
is either an eigenvalue, or a part of the continuous spectrum, see Ref. [ 13 , l3l| for a detailed analysis of a similar 
problem). This means that if Aoo G {to, to), then for a small change of parameters ao,&o the eigenvalue Aoo stays in 
the gap. We shall see one such case in Example 1111.61 


As already mentioned, m[o,nl is the [n/n + 1] Fade approximant of Wfo.oo) and it is an Afi function for n > 1. 
Consequently it can be represented in one of the forms (lll.5l) - (lll.7l) . As we have just done in Theorem 111. II for TO[o^oo), 
one can specify this representation: 

Proposition II.3 Each function m[o,n] (n = 1, 2,...) admits a representation 

m[o,n]{z) = - -—- ye (^ + dn+ [ , (n.l6) 

{z - Xn){z - Xn) V J t-z J 

with Xn G C“'", (i„ G M and pn a finitely supported measure. Furthermore, the parameters Xn and dn and the measure 
fin o,re uniquely determined and Xn is the unique eigenvalue of nonpositive type of -ff[o,ri] - 

The details of the proof of uniqueness can be found e.g. in [13, [13, IHi- The uniqueness in both Theorem III. II and 
Proposition III.3I guaranties that A„ (n = 1, 2,. ..) and Aoo are properly defined. In the literature they are called the 
generalized poles of nonpositive type of the corresponding Mi function, see [s^j- Using the classical result saying that 
TO[i_„] converges to TO[i_oo], see [^[^[S^j, one can prove -via eq. (111.3^ - the following convergence result, cf. [13]: 
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Proposition II.4 The functions TO[o „] (n S Z_|_) converge to m[o oo) as n ^ oo, locally uniformly on C+ \ ([^ 1 ,^ 2 ] U 
{Aoo}). 

Further generalization to different types of 7y-selfadjoint Jacobi matrices can be found in (^ . As a consequence we 
have the following corollary (cf. (H): 


Corollary II.5 The pole A„ o/m[o,„] converges to the pole Aoo ofm^Q oo) as n ^ 00 . 

Proof. If A € C+ or is an isolated eigenvalue then this statement is a simple consequence of Proposition IIL4I and 
the Rouche theorem. 

If instead Aoo is real and is not an isolated eigenvalue then from Proposition III. 41 and the Rouche theorem we see 
that all the accumulation points of A„ lie in the spectrum of iJ[o,oo)j which is a compact set. Since both the functions 
m[o,„] and TO[o^oo) belong to Afi and are of the type (III.6() . we have 


m[o,n]{z) = 


I 


{z - Xn)iz - An) 


^rTTn(z), 


^[0,oo 


(z) = 


1 


Aoo)(-2' Aoo) 




where <pn and ip are Nevanlinna functions and r„ and r are rational functions. Suppose now that there is a sub¬ 
sequence such that A„^ —>■ Aq Aoo- As a consequence, ipn^ should also converge to a Nevanlinna function ipQ ^ ip 
which contradicts the uniqueness oi (p. □ 


III. CONVERGENCE RATES 

Now we are in a position to ask the principal question of this paper: 

What is the character of the convergence of An Aoo ? 

We mainly consider the situation when Aoo is simple eigenvalue located outside the support of the measure fi in (III. 91) : 
we show a theoretical bound on the convergence rate and test it on examples. 


A. Theoretical results 


In this section we consider the situation when Aoo is a simple pole of m. Note that if Aoo G C"*", then it is necessarily 
a simple pole, due to Theorem lII.il (if: moreover, there exists no G N such that A„ € C+ for n > no (see Theorem lA.il 
below). If instead Aoo is a simple real pole, we show that A„ is real for sufficiently large n. We begin with a technical 
result, needed to prove our main theorem. 


Proposition III.1 Let m = nijo^oo) satisfy the (equivalent) conditions (i), (ii), (Hi) of Theorem If Aoo G 

C\ [^ 3 ,^ 4 ] is a simple pole of m, then 


Aoo An. — 


bl (mq 

,00) (-^oo ) ^[l,n](Aoo)) 


1-61 


0 [1.00) 


(Aoo) 


where an is such that 


sup,j,gx \‘mii,oo){x) - mii^oo)ix)\ 


n —>■ 00 


for any disc X C C"*" containing Aoo- If, additionally, Aoo G M, then A„ G M for sufficiently large n. 


Proof. Let X be an open disc, such that Aoo G A C C+ \ [^3,^4]. Note that for sufficiently large n the functions 

1 


i(z) := - = z - ao - nfz), 

mio,n]{z) 


(III.1) 


as well as 


W[o,oo)(^:) 


= z-ao- blm[i^oo){z) 


moo{z) 


(III.2) 
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belong to C(X), the complex Banach space of continuous functions on X with the supremum norm. Indeed, for 
sufficiently large n the function m[i „](z) has no poles in X fl K. Also observe that converges to rriao in C(X), 
since converges to TO[i_oo)(-z) locally uniformly on C \ [t 3 ,t 4 ]. Consider the mapping 

F : C{X) X X B {m,x) i—>■ m{x) G C. 

As Aoo G C"*" is a simple pole of m[o_oo)j one has 

^oo(’^oo) “ ( l/^[0,oo)) (Aoo) ^ 0. 

Therefore, 

dF 

Aoo) = ^oo('^) 7^ 

and we can apply the implicit function theorem in Banach spaces to the mapping F (see e.g. [1^). As a result we 
obtain in a neighborhood U xY oi (tooo, Aoo) a differentiable function ^ : U ^ Y such that 

{{m,x) G U X Y : m{x) = 0} = {(m,^(m)) : m G U} = 0. 

We may take Y so small that Y C C \ [tg, ^ 4 ] and that TOoo has no other zeros in Y except Aoo- Note that for sufficiently 
large n one has rrin G U. Hence, on one hand we have that for sufficiently large n 


mniiimn)) = F{mn, ^{mn)) = 0, mn{x) 0, x gU \ {^(m„)} . 


On the other hand, A„ converges to Aoo and mn{\n) = 0. Consequently, A„ = ^(mn) for n large enough. 
Now note that 


dx , , ^(moo,x(moc))m m(Aoo) 

-(moo)™ = “- 


dm 


dx 


(^007^(^00)) ^oo('^oo) 


Furthermore, 


dx 


•^oo '^n — Q (^^oo)(^cxD ^n) 


moo (A 00 ) -mn{X 00 ) 


™'oo(-^oo) 

where a(hn)/ ||^n|lc(x) ^ with —?► 0. Set = m„ — moo and 


+ a(moo - mn), 


— rr(moo ^n) — Aoo Ayi “b 


^oo(Aoo) ^n(Aoo) 
m'ooiXoo) 


and note that the right-hand side of the above does not depend on the initial choice of the disc X. This finishes the 
proof of the hrst statement. 

Now let Aoo G R \ [^ 3 , ^ 4 ]- By the locally uniform convergence of m„ to moo and by the Rouche theorem there is a 
small disc Z with the center in Aoo, such that each function mn{z) has precisely one zero Zn in Z. As A„ converges 
to Aoo we must have A„ = for large n. Therefore, An G K, otherwise Xn G Z is another zero of mn in Z, which is a 
contradiction. □ 


Remark III. 2 We are able now to prove the main result of our paper. Theorem IIII.3I First, though, we would like 
to stress that there are two equivalent ways of seeing it according to the objects we consider: 

A first interpretation takes as its main object the tridiagonal matrix, presented here in a block form 


Hi 


[0,oo) 


/ Qq 
bo 
0 

V i 


— bn 


0 




H, 


[ 1 . 00 ) 


Aoo is then the (unique) eigenvalue of nonpositive type of H[q_oo), A„ is the unique eigenvalue of nonpositive type of 
the finite truncation H[o^„] of H[o^oo)j a-nd the spectrum of Hq 00 ) is contained, by assumption, in [^ 3 ,^ 4 ]. 
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A second interpretation considers instead a meromorphic function m{z) having the representations (III. 81) and (III. 91) 
and its [n — \/n] Fade approximants The point Aoo (A„) is then the unique pole of nonpositive type of 771 ( 2 ) 

(77i[o,„], respectively). 

In both settings Theorem IIII.3I gives the convergence rate of A„ to Aoo , in terms of the “distance” of Aoo from the 
interval [^3,^4]: the rate of convergence is at least exponential where the number q is such that Aoo lies on the 

ellipse with foci at < 3,^4 and sum of its semi-axes equal to (ti — t3)q/2. Consequently, as confirmed by our numerical 
tests below, the convergence rate gets worse the larger is the eccentricity of said ellipse, i.e.: the convergence slows 
down when Aoo is “close” to the interval [^3,^4]. 

Theorem III.3 Let Xoo, A„, ti be as in Theorem mj\ and Remark \III.2\ above. If Aoo S C \ [^3,^4] is a simple 
eigenvalue, then 

limsup |Aoo - Anl^/” < 

n—¥oo Q 


where q = g + \Jg^ — and 


|Aoo — ^4] -I- 1 A 00 
ti — ts 



(III.3) 


is the reciprocal of the eccentricity of the ellipse through Aoo with foci at ^ 3 ,^ 4 . If, additionally, Aoo G R, then A™ G R 
for sufficiently large n. 

Proof. For i? > 1 let Ln denote the closed set bounded by the ellipse with foci at t^, t^ and the sum of its semi-axes 
equal to (^4 — tz)R/2. Due to Theorem (2.6.2) in Ref. [1^, one has 


limsup sup |771[1_„](2) - 7T7[i_oo)(2;)r^" < 

ri-J-oo zGC\Lr ^ 


(III.4) 


Note that 


Aoo ^ Tfl ^ R<g+ - 1 - (ni.5) 

Take any R G {l,g + g'^ — 1) and a small disc X, such that Aoo G X C <C \ Ln. From Proposition IIII.ll we obtain 

that 


|Aoo Ati I ^ |^[l,n] (*^00) ^[1,00) ('^00) I 0 !.ji 

< (^2 sup |77l[i^„](2) - m[i,oo)(2;)|, 

where Ci,C '2 are constants, dependent on i/[o,oo) and X only. As a consequence, 

limsup |Aoo - A„|^/" < 

n—foo 

Letting R^ g + \Jg^ — \ finishes the proof. □ 

Note that Theorem IIII.3I cannot be easily generalized to the case when Aoo S [tsjt^] \ supp/i. For example, if the 
support of the measure pi consists of two disjoint intervals [^ 3 , ts] U [tg, ^ 4 ] the estimate (IIII.4|) . which was the key point 
in proving Theorem IIII. 31 still holds only outside the ellipse with foci at is, t 4 , the reason being that the union of the 
poles of the Fade approximants of J^^it — z)~^fi{dt) may be dense in [^ 5 ,^ 6 ]) see for instance Ref. (3^ . 


B. Examples 


In our examples we want to be able to choose the position of Aoo; it is therefore convenient to consider cases where 
it is possible to calculate it without resorting at first to truncated matrices. One way to is look for matrices FI[i,oo) 
such that the the corresponding functions ttijj^ oo) have a closed, analytic form. This will allow us to use (numerical) 
root finding methods to calculate Aoo solving the equation 


2-00-60 W[i,oo)(2) = 0. 


(III.6) 
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Remembering that 


^[l,oo) 


Jt, t-z’ 


this reduces to finding a suitable measure fj,(t) with finite support. 
The choice 


dfi = daa,p(t) 


X[-i.i]W 


— s)“(l + sYds 


(IIL7) 


(IIL8) 


where (t) is the characteristic function of the interval [—1,1], gives us the matrix i7[i,oo) corresponding to Jacobi 

polynomials with parameters a,/3. To construct iJ[i,oo) "we consider the Jacobi orthogonal polynomials, which form 
an orthogonal basis in T^((Ta„s) and the multiplication operator p ^ xp in L‘^{aa,p)- three term recurrence 

relation (4.5.1) and the normalization factors (4.3.3) in Ref. [ij provide the tridiagonal representation iJ[i,oo) of the 
multiplication operator. The spectral theorem for selfadjoint operators guarantees that (IIII.7I) with (IIII.8I) is satisfied. 

The choice 


d^l = X[- 2 . 2 ] (0 • ^ J (ni.9) 

instead gives the matrix i7[i,oo) with aj = 0, bj = 1, j = 1,2 ,... corresponding to the orthogonal polynomials 
associated with the Wigner semicircle measure. 

To calculate A„ we use Matlab [l^: first, we calculate all eigenvalues and eigenvectors of i7[o,n]; we then find A„ as 
the only eigenvalue of nonpositive type of H[o,n] ^ i-e. the only eigenvalue which is either in the upper half-plane or is 
real with the corresponding eigenvector x satisfying 


n 

+ <0. 
i=i 

A summary of the relevant results for all our examples can be found in Table HI Only graphs useful to our discussion 
are shown here; for the remaining cases quoted, pictures can be found in the supplementary files [3^ : the file names 
there refer to those given in Table HI 

Example III.4 As our first example we take a = ,8 = 0 in eq. (IIII.8I) so that 

^ ^ - log(-l - z)), (III.IO) 

where the branch of the logarithm is chosen in such way that the above function is a Nevanlinna function. In this case 
H[i,oc) corresponds, in the way described in the remarks above, to the Legendre polynomials. We now vary the only 
remaining free parameters Og and bg of i7[o,oo)- Note in any case, due to Remark III. 21 we have X^o ^ [^ 3 )^ 4 ] = [~li !]• 
Let us start with ag = 0.5, 6o = 0.05. It is immediately evident from Figure[T]that the points A„ arrange themselves 
on three branches which spiral into Aoo. The point A„ jumps in a regular fashion from one branch to another: all 
A„’s with n mod 3 = const fall on the same branch. The two branches starting on the real axis leave it at n = 45 
and n = 190 respectively, as can be seen from the plot of the imaginary part of A„ in Figure [21 the plots of the 
corresponding real parts have each a cusp at the same time, due to the inversion of the direction of motion of A„. The 
plot in Figure [3] shows exponential convergence of A„ to Aoo setting in soon after all three branches leave the real line. 
The black reference line with slope —2 log(q) represents the bound from Theorem lIIL31 the intercept is chosen so that 
the line is superimposed on the numerical data. In this example (as well as in subsequent examples with measures 
having no gaps) we see that the estimate of the convergence rate in Theorem IIIL3I is sharp and consequently can not 
be improved in general. On the other hand, the estimate is not always sharp, cf. Example IIII. 61 

We have already mentioned that the points A„ arrange themselves regularly on branches. This behavior is common 
to most of our examples, except the cases when Aoo is on the real axis, where because of Proposition IIII. II there is 
one branch only. Since the branches appear to approach Aoo isotropically, as can be e.g. seen in the zoomed picture 
of Figure [U this suggests a convenient way of calculating the numerical value of Aoo as the mean 

N 

= k-^ ^ A„, 

n=A^ —fc+1 


(iii.ii) 
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where k denotes the number of branches and N the last n for which A„ is calculated. 

The value of Aoo obtained taking the average of the three last points (one for each branch) equals 0.4999 + 0.00391, 
which agrees well with the value 0.498631 + 0.003913971 obtained solving with Mathematica eq. (IIII.6I) . 

We conclude noting that -here and in all our other examples- the real axis forms a barrier for A„: it never crosses 
it and branches which touch it get stuck in it; this is clearly visible in the movie Legendre_spirals_a_0.5.avi that can 
be found among the supplementary files to this paper . This behavior is related to the symmetry of the spectrum 
with respect to the real axis. 

If we now keep oq = 0.5 constant and vary 6o to assume the values 6o = 0.5,0.1, and 0.01, in all cases the A„’s 
arrange themselves over three spiraling branches and the value Aoo obtained from the average of the last three values 
of Xn agrees with the one calculated solving eq. (IIII.6|) with Mathematica to the last digit shown in Table m In all 
cases convergence is exponential and the slope of logarithmic plots similar to that in Figure [3] are in agreement with 
the value —21og((3') from Theorem IIII. 31 fsee TableH]). There are some case to case differences, but they do not affect 
the picture given above: for 6o = 0-5 the three branches are not clearly visible, due to the very fast convergence of 
An; for 6o = 0.1 only one branch spends some time on the real axis, up to n = 46; and for b = 0.01 even after 1000 
iterations one of the three branches has not yet left the real axis (for figures see the supplementary material [4l|). 

Our last example with measure eq. (IIII. 101) is a case when Aoo G K \ [^3,0]: if we take oq = 1.001 and bo = 0.001 
we get Aoo — 1.001 (Mathematica has problems solving eq. pil.6|) in this case). The values of An are all real; we 
therefore have a single branch. Convergence is again exponential (for figures see the supplementary material |42j|b 
It is instructive to compare this case to the case ag = 0.5, b = 0.05: the distance of Aoo from the interval [—1,1] is 
of the same order, but the convergence is much faster in the present case. This is due to the different eccentricities 
of the ellipses from Theorem IIII. 31 in the present case the eccentricity is smaller, and therefore q is larger and in 
consequence the convergence rate is better. 

Finally, in Figure U we summarize the convergence behavior when submatrix i7[i,oo) corresponds to measure eq. 
(IIII. 101) : we vary the parameters ag € (—1,1) and bg € (0.05, 0.5) and for each pair (ag, bg) we compute the value no 
where the last A„ branch leaves the real axis. Each pair (ag, bg) determines a single point Aoo in the upper half-plane, 
which we plot color coded according to ng. The figure appears to have a fractal character for which we do not yet have 
an explanation but which we suspect to be related to the way the number of A„ branches varies with varying ag. This 
can be seen in the movie Legendre_spirals_b_0.01.avi where we keep bg = 0.01 and vary ag. We instead observe no 
change in the number of branches when varying bg at constant oq (see e.g. the movie Legendre_spirals_a_0.5.avi [^b 

As we have already mentioned, branches are a common occurrence, not limited to the example just given. We give 
here a couple more examples. 

Example III. 5 We now take ff[i^oc) corresponding to the orthogonal polynomials associated with the Wigner semi¬ 
circle measure, i.e. the measure given by eq. (IIII. 91) and we choose ag = 0.5, bg = 0.1 as the remaining parameters 
for i7[o,oo)- Note due to Remark (11.21 we again have Aoo ^ [^ 3 ,^ 4 ] = [—2,2]. 

In FigureOit is possible to see that the A„’s form twelve branches. Knowing this, we can use the recipe given above 
to calculate Aoo as the mean eq. (IIII. Ill) of the last twelve values of A„; the result agrees to the last digit shown in 
Table U with the value obtained solving eq. (IIII.6|) with Mathematica. 

Figure|6]shows exponential convergence of A„ to A^ with the rate predicted by Theorem lIII.3l other plots concerning 
this example can be found as supplementary files [43|. 

The video Wigner_spirals_b_0.01.avi in [s^ shows the evolution of the A„ branches under the change of the parameter 
ag. As in Example IIII. 41 it is evident that here too the number of branches changes with oq. Looking at the first 
few frames of the movie it is also evident that there are branches that start off the real axis, hit it and -instead of 
continuing into C“- get trapped in it moving horizontally for a number of n, and then leave it when the spiral reenters 
C+. 

Example III. 6 Here we present an example of a different nature. The matrix i7[i,oo) is constructed in such way, 
that its spectrum is a totally disconnected, Cantor-like set, see for details; the other parameters are ag = 0.5 and 
60 = 0.1. In this particular case we could count 27 branches of A„; such a high number is hardly visible when plotting 
An in the complex plane but can be seen in the plot of log |A„ — Aoo| in EigurejT] the A„’s form regular clusters of 27 
points. 

What is particularly noteworthy is that -contrary to the previous examples- the convergence rate is faster than the 
one predicted by Theorem IIII. 31 as can be seen comparing the theoretical bound (black line) with the numerical points 
in FigurejT] This is probably due to the fact that the spectrum of i7[i 00 ) contains gaps. Other pictures corresponding 
to this example (called Cantor_a_0.5_b_0.1_>oo<.eps), as well as the video Cantor_spirals_b_0.2.avi showing the evolution 
of branches under the change of the parameter ag, can be found in the supplementary files [ 1 ^: again the number 
of branches changes with oq; moreover intervals in ag where both Aoo and all the A„ become real are evident and 
correspond to the gaps in i7[i,oo)j see Remark (11.21 
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Although Theorem nil. 31 proves the exponential rate of convergence of A„ to Aoo, it does not say when this conver¬ 
gence starts manifesting itself: at least in theory we could have |Aoo — A„| > for some n. Our numerical tests 
indicate that the convergence is somehow “better” than the above bound, as the asymptotic behavior is K ■ < 7 “^” with 
K < 1. Even the curious initial behavior connected with the real axis we observed in several cases above, does not 
bring |Aoo — A„| to exceed q~^^. On the other hand it would be convenient to have an estimate of no such that for 
n > no the point An is sure to be outside the support [t3,t4] of the measure fi: if one or more of the An’s branches is 
still on [to, 14 ], our justification for estimating Aoo by the mean eq. (IIII.llI) is compromised. We give an upper bound 
for no in Appendix [Al 


C. Aoo € R is embedded in the spectrum of the representing measure 


We shall limit our investigation of the case when Aoo G [^ 3 ,^ 4 ] to an example where Aoo = ts, to show that 
convergence in this case is in general worse than exponential. To build our example, we start by recalling Remark 
III.21 The only example known to us of classical orthogonal polynomials satisfying condition (III.lSp with Aoo S [t 3 ,t 4 ] 
are the Jacobi polynomials with parameters a > 2 ov f3 > 2 and with Aoo = to or Aoo = t 4 , respectively. 

Example III.7 We take JJ[i,oo) such that 


^[1,00) 





i.e. we take a = 2,/3 = 1 in eq. (IIII. 8 I) . We then choose ao = —5/3, bo = \/2/3, so that Aoo = —1- This can be 
seen e.g. by analyzing the matrix iJ[o,oo) itself: it results from conditions (i) and (ii) that —1 is an algebraically 
simple eigenvalue of iJ[o,oo) with the corresponding eigenvector x satisfying —|a;op + = 0 ; i- 6 - "^00 = —1 is 

the unique eigenvalue of nonpositive type of 

Furthermore, due to Theorem 2.2 (pis) of [2lf, —1 is a singular critical, algebraically simple eigenvalue (see [2lj for 
a classification of eigenvalues of nonpositive type). 

The log-log plot of |A„ — A| Vs. n in Figure [ 8 ] shows clearly that convergence in this case is only polynomial: 
|Ara — A| ~ as can be seen comparing the numerical data with the black reference line whose slope is —2. The 
plot of the real and imaginary part of A„ can be found in [s^ as the file Jacobil2crit_Relm.eps. 


IV. CONCLUDING REMARKS 

As already stated above, the main question we tried to answer here is: “Suppose we are only able to compute 
the spectrum of finite truncations of a pseudo-Hermitian tridiagonal matrix iJ[o,oo); can we say something about the 
spectrum of the full operator iJ[o,oo)? Most interestingly, can we predict if the spectrum of iJ[o,oo) is real or not?” 

Suppose we have performed N iterations, increasing step by step the size of the truncated matrix JJ[o,n]; one of the 
next four cases applies. 

• An is real for all n = 1,...,V, is either the minumum or the maximum of the spectrum of J?[o,n] arid is 
separated from the other eigenvalues. Then the limit point Aoo is also a real single eigenvalue, separated from 
other eigenvalues of iJ[o,oo)- 

• An is complex for all n larger than a no < N. Then the limit point Aoo is also a complex single eigenvalue. Aoo 
itself can be evaluated by finding the number of An branches and then using eq. (IIII.llI) . 

• An oscillates between the real line and the complex plane up to n = V, a common occurrence when Aoo is very 
close to the support of the measure jjL. Then the situation is in principle unclear: the limit eigenvalue Aoo might 
be a complex point, a real critical point, or a real point in a relatively small gap of the spectrum. Still, if the 
An branches can be found and not to many of them are still trapped on K for n ~ V, it is still possible to 
give a numerical evaluation A^^ of Aoo by a careful use of eq. (IIII.lip . If the plot of log I An — A^^l is then 
approximately a straight line, this is another indication for Aoo being a simple eigenvalue out of the support of 

• The sequence An converges to a point Aoo £ K, but the convergence is not exponential, which again can be 
seen by the study of the plot of log|A„ — Aoo|- Then Aoo is a critical point on the real line, embedded in the 
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support of In view of the second equation of (III.151) . this case seems to be non-generic, i.e. a small change 
of the entries of the matrix will lead to a different case. However, we were able to clearly observe this case in a 
numerical simulation, which in our opinion is an argument for considering this possibility as well. 

While, for sake of simplicity, we restricted ourselves to the case of matrices with a single eigenvalue of nonpositive 
type, we believe that the results derived above can be generalized to a wider class of operators, at least to those 
considered in Ref. [^ . 

In the context of random matrices or Nevanlinna functions, our research can be viewed as concerning the 
problem of predicting whether or not Aoo G [ts, ^ 4 ] by calculating a finite number of Fade approximants. A connection 
can also be found to Fade approximation of the Z-transform, considered in , where the real line is replaced by the 
unit circle. 
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Appendix A: Estimate of the maximum number of real A^, when Aoo G 
Theorem A.l Let m = TO[o,oo) 0 / the forms (|II.8I) and (III.91) with Aoo € C’*'. Then for any e G (0,1) for 

naxjggtp |m[o,o 


n > := 


f |m[o,oo)(z)l 

I (g-l)Ut 4 -* 3 ) 

log (l + ey^l- 


(A.l) 


has 


A„ G C \ 


where g is again given by eq. (IHI.3I) . i?(e) = g + ey/g^ — 1 and again denotes the closed set bounded by the 

ellipse with foci at ^ 3,^4 and the sum of its semi-axes equal to (fA — t^)R{£)/2. In particular, if n > for some 
£ G (0,1) then A„ ^ [^3,0]- 


Proof. Consider the functions 


m„{z) := - iriooiz) = 


We now recall the estimate given in as formula (6.10): for every n G and for every i5 G (1,R) one has 

8,5 


Sn = S (I, _ tMS - 1)2 Vfl 


(A.2) 


Applying equation (IA.2I) with d = g, R = R{£) = g + £\/ <7^ — 1 we get after elementary transformations of (jA.lIl 
that 


|m„(z) - m^{z)\ 


= &ol™[l.n](^) - W[i_oo)(2)l 
^ ( 9 

{tA-U){g-lf yiA^e^g2_i 


< 


1 

max^^zSLa |m[o,oo)(^)| 


min |moo(2;)|. 
z£SLr 
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Hence, by the the Rouche theorem, m„ and mao have the same number of zeros in C\ L/j. However, TOqo has precisely 
two zeros in C \ L/j, namely Aoo and Xao- As A„ is the only a zero of in the upper half-plane, we get A„ € C \ L/j. 

□ 

If we now apply Theorem lA.ll to Examole lllLdl we get A^o .5 = 52 for bo = 0.5, A^o .5 = 2155 for bo = 0.1, TVq.s = 10917 
for bo = 0.05, and TVg.s ~ 4 • 10® for bo = 0.01. Comparing with Example IHI.4[ where no — 1,46,190, and uq >> 1000 
respectively, we see that the estimate N 0.5 is a far from tight upper bound whose ratio to the numeric value no appears 
to grow with decreasing bo¬ 
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oo 

ho 

maxn 

Acxd 

no 

q 

1 Aoo — An 1 

Legendre 

0.5 

0.5 

50 

0.4045 -h 0.3064i 

1 

1.3855 

~ 


0.5 

0.1 

200 

0.4946 -h 0.0155i 

46 

1.0180 



0.5 

0.05 

1000 

0.4986 -h 0.0039i 

190 

1.0045 

~ 4.6-ig-2" 


0.5 

0.01 

1000 

0.4999 -h 0.0002i 

> 1000 

1.0002 

~ 6.5"^g"^" 


1.001 

0.001 

200 

1.0010 

— 

1.0456 

~ 12.5"^g“^’" 

Wigner 

0.5 

0.1 

1000 

0.4975 -h 0.0096i 

175 

1.0050 


Cantor 

0.5 

0.1 

400 

0.5074 -h 0.0096i 

81 

1.0044 


Jacobi(l,2) 

-5/3 

72/3 

300 

-1 

— 

1 

~n-2 


TABLE I: Parameters and numerical results for Examples IIII.4I [III.51 HILHI and IIII.7I If Aoo ^ R then no denotes the maximal 
n for which An £ R, ? was computed according to Theorem IIII.3I 
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FIG. 1: The only eigenvalue An of nonpositive type of the matrix lF[o,n] given in eq. (11.21) with ao = 0.5, bo = 0.05 and 
being the Jacobi matrix corresponding to the Legendre polynomials. The points, each corresponding to a different n, 
are plotted on the complex plane color coded according to n in the upper picture, and according to n mod 3 (0-blue, 1-cyan, 
2-yellow) in the lower picture, where we show a zoomed detail around Aoo. The three branches of An mentioned in the text are 
clearly visible. See Example IIII. 41 













17 



FIG. 2: Real and imaginary part of the eigenvalue An from Figure [T] Vs. n. The points color coded according to n mod 3. 
The cusps in the real parts of the yellow and cyan branches correspond to the change of direction when a branch leaves the 
real line, see Figure [T] See Example IIII. 41 



FIG. 3: Same parameters as in Figure [T] The black line represents the theoretical bound given in Theorem IIII. 31 shifted down 
so as to be superimposed on the numerical data. The logarithmic scale on the vertical axis makes evident the exponential 
convergence of An to A OO • See Example IIII. 41 
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FIG. 4: The value no where the last An branch leaves the real axis as a function of Aoo- l^[i,n] is the Jacobi matrix corresponding 
to the Legendre polynomials and ao and bo vary to get the different values of Aoo- See Example IIII. 41 



Re(A) 


FIG. 5: The eigenvalues An for IF[i,n] associated with the Wigner measure and parameters ao — 0.5, bo = 0.1, plotted in the 
complex plane, color coded according to n mod 12 where twelve is the number of An branches. See Example IIII. 51 
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FIG. 6: Same parameters and color code as in Figure [5] The black line represents the theoretical bound given in Theorem 
nil. 31 shifted down so as to be superimposed on the numerical data. The logarithmic scale on the vertical axis makes evident 
the exponential convergence of A„ to Aoo. See Example IIII. 51 



FIG. 7: Exponential convergence of A^ to Aoo. Here H[i^ao) is associated with a measure on a Cantor-like set and ao = 0.5, 
bo = 0.1. The black line represents the theoretical bound given in Theorem IIII. 31 Color coding according to n mod 27 
evidences a repeated pattern of 27 points. The observed convergence is again exponential but much faster than the bound 
given by Theorem IIII. 31 the observed g is approximately 1.1147, much larger than the theoretical one 1.0044. See Example 

ULi 
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FIG. 8: Log-log plot showing 0{n convergence of An to Acx^ = —1 in a case when Aoo is embedded in the spectrum of L/p 
See Example nil. 71 






